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The Tevatron has inspired new interest in the subject of magnetic monopoles. First 
there was the 1998 DO limit on the virtual production of monopoles, based on the theory 
of Ginzberg and collaborators. In 2000 the first results from an experiment (Fermilab 
E882) searching for real magnetically charged particles bound to elements from the CDF 
and DO detectors were reported. This also required new developments in theory. The 
status of the experimental limits on monopole masses will be discussed, as well as the 
limitation of the theory of magnetic charge at present. 

1. Maxwell's Equations 

The most obvious virtue of introducing magnetic charge is the symmetry thereby 
imparted to Maxwell's equations, 

V-E = 47rpe, V-B = 47rp„, (la) 

„ Id Att. „ 1 9 „ 47r . 

VxB = --E+— je, -V X E= --B + — j„. (lb) 
c ot c c ot c 

These equations are invariant under a global duality transformation. If £ denotes 
any electric quantity, such as E, Pf,, or je, while M denotes any magnetic quantity, 
such as B, pm, or j,„, the dual Maxwell equations are invariant under 

£ ^ Scose + Msm0, M ^ McosO - £ sin 6*, (2) 

where is a constant. 

J. J. ThomsonS (1904) observed the remarkable fact that a static system of 
an electric (e) and a magnetic {g) charge separated by a distance R possesses an 
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Fig. 1. The relative motion of two dyons is confined to the surface of a cone about the direction 
of the angular momentum. 



angular momentum, 

E X B 



J 



j {dr) rxG = J (dr) r x 



4ttc 



er q(r — R) 

-TT X 



(r-R) 



^R, (3) 



which follows from symmetry (the integral can only supply a numerical factor, which 
turns out to be At:). The quantization of charge follows by applying semiclassical 
quantization of angular momentum: 

J-R= — =n-, or eg^-hc, n = 0, ±1, ±2, . . . . (4) 
c 2 2 



2. Classical Scattering 

Actually, earlier in 1896, Poincarei investigated the motion of an electron in the 
presence of a magnetic pole. Let's generalize to two dyons (a term coined by 
Schwinger in 1969) with charges ei, gi, and 62, 52, respectively. There are two 
charge combinations q — 6162 + 5i(?2, « = — .£i22-i-£Mi , Then the classical equation 
of relative motion is {fi is the reduced mass and v is the relative velocity) 

/^lir^V-^vx-^. (5) 
The constants of the motion are the energy and the angular momentum, 

E = \^iv^ + - , J = r X + Kf . (6) 

2 r 

Note that Thomson's angular momentum is prefigured here. 

Because J • f = k, the motion is confined to a cone, as shown in Fig. |]. Here the 
angle of the cone is given by 

X ^ 

cot - = — , 1 = fivob, (7) 
2 \k\ 
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where vq is the relative speed at infinity, and b is the impact parameter. The 
scattering angle 9 is given by 



X 

cos - = COS — 

2 2 



tan 



e/2 
cosx/2 

1 f k|"o 



where 



q cot f y , 



q>0, 



TT — tan 



-1 / \k\vo 



cot , gr < 0. 



(8a) 
(8b) 



The impact parameter b{9) is a multiple- valued function of 9. The differential cross 
section is therefore 



da 



bdb 



d{cos9) 



E 



1 



2^vo J sin 



^4 X 



sinxdx 



sin 9 d9 



(9) 



Representative results are given in Ref. 0. 

The cross section becomes infinite in two circumstances; first, when 

sin6l = (sin X 7^0), 9^?:, (10) 

we have what is called a glory. For monopole-electron scattering this occurs for 



^ = 1.047, 1.318, 1.403, .... 

The other case in which the cross section diverges is when 

d9 



0. 



(11) 



(12) 



This is called a rainbow. For monopole-electron scattering this occurs at 

9r = 140.1°, 156.7°, 163.5°, .... (13) 
For small scattering angles we have the generalization of the Rutherford formula 



da 

dn 



1 



6132 - 62^1 



6162 + gig2 



1 



{9/2Y 



61 <1. (14) 



3. Quantum Theory 

DiracI showed in 1931 that quantum mechanics was consistent with the existence of 
magnetic monopoles provided the quantization condition (^ holds, which explains 
the quantization of electric charge. This was generalized by Schwinger to dyons: 

71 

6l52 - 6231 = -flC. (15) 

(Schwinger sometimes argued that n was an even integer, or even 4 times an integer.) 
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One can see where this conies from by considering quantum mechanical scatter- 
ing. To define the Hamiltonian, one must introduce a vector potential, which must 
be singular because B 7^ V x A. For example, a potential singular along the entire 
line n is 

A(r)^-glf - U-^cot00 if n = z, (16) 



r 2 \r — h - V r + n- r/ r 

which corresponds to the desired magnetic field from a magnetic monopole, B(r) = 
g^. Invariance of the theory (wavefunctions must be single- valued) under string ro- 
tations implies the charge quantization condition. This is a nonperturbative state- 
ment. 

Yang offered another approach, which is fundamentally equivalent^ He insisted 
that there be no singularities, but rather different potentials in different regions: 

v4;^ = (1 -cos6') = -tan^, 6* < tt, (17a) 

^ rsmt^ r 2 

A\ = (1 + COS0) = -^cot^, ^>0. (17b) 

^ rsmt^ r 2 

These correspond to the same magnetic field, so must differ by a gradient: 

Requiring now that e^'^'^ be single valued leads to the quantization condition (^. 

There is also a intrinsic spin formulation, pioneered by GoldhaberJj The energy 
(^ differs by a gauge transformation from 

1 ( 2 . J'-(J^f)' ^ , eie2+gi.g2 

where 

S X r 

J = rxp + S, ^v = p + — (20) 

The quantization condition appears as S • f = m'. This was elaborated long agoE 
The nonrelativistic Hamiltonian for a system of two interacting dyons is 

TiV^s- 2m' cos 61 1 a m'^ „\ q 

W = -— + ^-TT- ^cot^ei -f-i, 21 

2[i \ r2 sin^ Q i d<^ J r ^ ^ 



where m' — ~{eig2 — e2gi)/hc. The wavefunction separates: V'(r) = i?(r)6(0)e""'^, 
where 



(P 2d 2iiq 7(7 + 1) - m'2 > 



e-j(j + i)e. (22b) 



1 


d 




_ sin I 


)d9 





sin^ ( 
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The solution to the 9 equation is the rotation matrix element: {x — cos 9) 



where p^™'"' are the Jacobi polynomials, or "monopole harmonics." This forces m' 



j 

to be an integer. The radial solutions are confluent hypergeometric functions, 

Rkjir) ^ e-'''''{kr)^F{L + l-iT],2L + 2,2ikr), (24a) 



2 



We solve the Schrodinger equation such that a distorted plane wave is incident, 
i/jin = exp i [k • r + r/ ln{kr — k • r)] . (25) 
Then the outgoing wave has the form 

V-out ~ }.eiikr-r,ln2kr)^^m'^' ji^g^^ ^26) 

r 

so up to an unobservable phase, the scattering amplitude is (here 9 is the scattering 
angle) 

oo 

2^kfi9) = J2 (2J + - 9)e-^^^^-'''^\ (27) 

j=\m'\ 

in terms of the Coulomb phase shift, Sl = argr(L + 1 + irj). Note that the integer 
quantization of m' results from the use of an infinite ("symmetric") string; an 
unsymmetric string allows m' = integer + i . 

One can show that reorienting the string direction gives rise to an unobservable 
phase.Ll Note that this result is completely general: the incident wave makes an 
arbitrary angle with respect to the string direction. Rotation of the string direction 
is a gauge transformation. 

Notice that small angle scattering is still given by the Rutherford formula: 



da ( m' 



dfl \2k J sm^9/2 



« 1, (28) 



for electron-monopole scattering. The classical result is good roughly up to the first 
classical rainbow. In gen£ral, one must proceed numerically. Various remarkable 
results are shown in Ref. cl. 

We can also include the effect of a magnetic dipole moment interaction, by 
adding a term to the Hamiltonian, 

-|^7'T-H, n = g^. (29) 
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Relativistic Cross Sections 

Kazama, Yang, and Goldhaber 




90.0 110.0 130.0 150.0 170.0 

9 



Fig. 2. Relativistic scattering of a spin-1/2 electron by a heavy monopole. 



For small scattering angles, the spin-flip and spin-nonflip cross sections are (for 
7= 1, 1) 



da 

dh 



m 
2k 



sin^ 9/2 


da 






sin* e/2' 


dn 


NF 





■9/2 



(30) 



Note that the spin-flip amplitude always vanishes in the backward direction; the 
spin-nonflip amplitude also vanishes there for conditions almost pertaining to an 
electron: for m' > 0, and 7 = 1, ^n^"") ~ Various results are shown in Ref. i. 

All of this work was done many years ago in Ref. i, which just goes to show that 
"goocLwork ages more slowly than its creators." There was of course much earlier 
work cIeI a relativistic calculation of the scattering of a spin-l/2_Dirac particle by 
a heavy monopole was given by Kazama, Yang, and Goldhaber.113 They also gave 
helicity-flip and helicity-nonflip cross sections which are shown in Fig. ^. (Note for 
9 = n, helicity nonflip corresponds to spin flip, and vice versa.) 



4. Quantum Field Theory 

The quantum field theory of magnetic charge has been developed by many people, 
notably SchwingeiE] and ZwanzigerEl A recent formulation suitable for eikonal 
calculations is given in Ref. [1. Formal Lorentz invariance is demonstrated provided 
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the quantization condition holds (rationaUzed units): 

Baffb - ebflg _ 1 f , unsymmetric 1 ^ (-^i) 

4tt \ n , symmetric [ ' 

( "Symmetric" and "unsymmetric" refer to the presence or absence of dual symmetry 
in the solutions of Maxwell's equations.) 

The electric and magnetic currents are the sources of the field strength and its 
dual: 

d-'F^, = and 8" *F^, = % , (32) 

where 



1 

which imply the dual conservation of electric and magnetic currents, jf^ and 
respectively, 



d^f = and 9^ *f = . (34) 

The first-order form of the action describing the interaction of a spin-1/2 electron 
■0 and a spin-1/2 monopole x is 

W = J{dx)(^ - ^F^'•'{x) {d^A, (x) - d,A^ (x)) + ^F^,{x)F^''' {x) 

+ '4!{x) {ijd + ejA{x) - m^) jp{x) + x(x) {i^d + gjB{x) - m^) x(a;)| , (35) 

where and F^^i, are independent field variables and 

Bf.{x) = ~J {dy)r {x - y) *F,,{y) , (36) 

Here ffj.{x) is the Dirac string function which satisfies the differential equation 

d,r{x) - Six) , (37) 

a formal symmetric solution of which is given by 

nx)^n^in-d)-'6ix), r(a;) = -^(-x), (38) 

where is an arbitrary vector. A corresponding dual form in terms of independent 
variables and *F^i, can be immediately written down. AlthougiLfrom these 
actions a complete path integral version of dual QED can be givenjl3 all that is 
needed for our purposes here is the relativistic interaction between spinor electric 
and magnetic currents = e'tp^^^'ip and = gxitiX- 

WU, *j) = I {dx){dx'){dx") Tix)e,,,rdT {x - x') D+ [x' - x")f {x") , (39) 
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where the photon propagator is denoted by D^{x ~ x'). 

The modern path integral reformulation of the dual quantum electrodynamics 
of electric and magnetic charges were used in Ref. 113 to rederive and generalize 
the eikonal results of Urrutiao for high-energy, low momentum-transfer scattering 
between an electron and a monopoie. A simplified version of the argument appears 
in Ref. E3. In terms of the quantization condition the scattering amplitude, for 
an arbitrary direction of the incident particle, turns out to be 

/(q) = -^e-2»^ (40) 

q 

where q is the momentum transfer. Squaring this and putting in the kinematical 
factors we obtain Urrutia's resultEJ 

dt At: t^' ^ ' ^ ' 

which is exactly the same as the nonrelativistic, small angle result found in Eq. (|2S 
This calculation, however, points the way toward a proper relativistic treatment, 
and will be extended elsewhere to the crossed process, the production of monopole- 
antimonopole pairs through quark-antiquark annihilation. 

5. Previous Searches for Magnetic Monopoles 

In the context of "more unified" non-Abelian theories, classical composite monopoie 
solutions were discovered. The mass of these monopoles would be of the order of 
the relevant gauge-symmetry breaking scale, which for grand unified theories is of 
order 10^^ GeV or higher. But there are models where the electroweak symmetry 
breaking can give rise to monopoles of mass ~ 10 TcV. Even the latter are not yet 
accessible to accelerator experiments, so limits on heavy monopoles depend either 
on cosmological considerationsjl3 or detection ol cosmologically produced (relic) 
monopoles impinging upon the earth or moon.llZI However, a priori, there is no 
reason that Dirac/Schwinger monopoles or dyons of arbitrary mass might not exist: 
In this respect, it is important to set limits below the 1 TeV scale. 

At the University of Oklahoma we are carrying out a direct search for monopoles 
produced at the Tevatron, which we will describe in the next section. But indirect 
searches have been proposed and carried out as well. De RujulaES proposed looking 
at the three-photon decay of the Z boson, where the process proceeds through a 
virtual monopoie loop. If we use his formula for the branching ratio for the Z — s- 87 
process, compared to the current experimental upper limitlij for the branching ratio 
of 10~^, we can rule out monopoie masses lower than about 400 GeVjrather than 
the 600 GeV quoted by De Rujula. Similarly, Ginzburg and PanfilE2l and more 
recently Ginzburg and Schilleitil considered the production of two photons with 
high transverse momenta by the collision of two photons produced either from e^e" 
or quark- (anti-) quark collisions. Again the final photons are produced through a 
virtual monopoie loop. Based on this theoretical scheme, an experimental limit 
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has appeared by the DO collaboration,Ej which sets the foUowing bounds on the 
monopole mass M: 

r 610 GeV for 5 = 
> { 870 GeV for 5* = 1/2 , (42) 
1580 GeV for = 1 



n 



where S is the spin of the monopole. It is worth noting that a lower massJimit of 
120 GeV for a Dirac monopole has been set by Graf, Schafer, and Greiner,E3 based 
on the monopole contribution to the vacuum polarization correction to the muon 
anomalous magnetic moment. (Actually, we believe that the correct limit, obtained 
from the well-known textbook formula for the g-factor correction due to a massive 
Dirac partkile is 60 GeV.) 

In Ref. Ila we have criticized these limits on theoretical grounds. They are based 
on a naive application of duality, in which the quantization plays no role. Thus gauge 
invariance is not demonstrated. The Euler-Heisenberg Lagrangian is used outside its 
range of validity for hard photon processes. That the Eulcr-Heisenberg Lagrangian 
is not an effective Lagrangian in the sensa^ol capturing radiative corrections is 
demonstrated by the disparate work of Refs. Ej'tJ. Moreover, the substitution e — > 
or 

137 

a — > Qfg = ' ^ = li 2, 3, . . . , (43) 

is made, which implies the manifest inconsistency of perturbation theory (which is 
already precluded by the nonperturbative quantization condition). The expansion 
parameter is Ofg, which is huge. Instead of radiative corrections being of the or- 
der of a for the electron- loop process, these corrections will be of order ag, which 
implies an uncontrollable sequence of corrections. For example, the internal radia- 
tive corrections to the four-photon box diagram have been computed by Ritufl 
and by Renter, Schmidt, and Schubert0 in QED. In the 0{a^) term in the ex- 
pansion of the EH Lagrangian, the coefficients of the (F^)^ and the (F*F)'^ terms 
are multiphed by (l f f -h 0{a^)) and (l + -I- 0{a'^)) respectively. These 

corxfictions become meaningless when we replace a ag. Moreover, it is easy to 
seeta that unitarity is violated by the formulas used in the DO analysis unless the 
monopole masses are above 1 TeV, so the limits quoted are meaningless. 

6. Oklahoma Experiment: Fermilab E882 

The best prior experimental limit oil the direct accelerator production of magnetic 
monopoles is that of Bertani et alS: a < 2x lO^^^cm^ for M < 850 GeV. The 
fundamental mechanism is supposed to be a Drell-Yan process, p + p m + rfi + X, 
where the cross section is given by 



= (68.5n)^^38™! f ^ ^ 
dM 9s J xi \ sx\ J 



Here is the invariant mass of the monopole-antimonopole pair, and we have 
included a factor of (3^ to reflect (1) phase space and (2) the velocity suppression of 
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the magnetic coupling. Note that we are unable to calculate the elementary process 
QQ ^ 1* ^ mfh perturbatively, so we must use nonperturbative estimates. 

Any monopole produced at Fermilab loses energy as it passes throueh the de- 
tector by ionization (computed using the energy loss formula of Ahlcnca) and is 
trapped in the detector elements with 100% probability due to interaction with 
the magnetic moments of the nuclei. The experiment consists of running sam- 
ples obtained from the old DO and CDF detectors through a superconducting in- 
duction detector. A schematic of our induction detector is shown on the web at 



tittp : //www. nhn. ou. edu//o7Egrk/apparatus .pdf . We are able to set much better 
limits than Bertani et al. because the integrated luminosity delivered to DO is 10^ 
larger than that of the previous 1990 experiment: J C = 172 ± 8pb~^. 
If (7 = 6162 + 5i(?2 < 0, TYnr in Eq. (21) gives binding of dyons 



+ - + ((.7 + 1/2)2- 



(45) 



Monopoles will not bind this way — a magnetic moment coupling as in Eq. (29) is 



required, in terms the gyrornagnetic ratio 7 = 1 + k = |. (7 = 1 or g = 2 is the 
"normal" value.) The theoryt£l is somewhat complicated and most inconclusive, as 
Table |l| shows. 



Table 1. Weakly bound states of nuclei to a magnetic monopole. The angular momentum quantum 
number J of the lowest bound state is indicated. In Notes, NR means nonrelativistic and R 
relativistic calculations; he indicates an additional hard core interaction is assumed, while FF 
signifies use of a form factor. IM means induced magnetization, an additional interaction employed 
for the relativistic spin-1 calculation. We use n = 1 except for the deuteron, where n = 2 is required 
for binding. 



Nucleus 



Spin 



7 



J 



Eb 



Notes 



Ref 




n 



2 He 
?IA1 
flAl 
il^Cd 



-1.91 
2.79 



0.857 
-2.13 
3.63 
3.63 
-0.62 



I - 1 
I 



I - 



(n = 2) 

_ 3 
2 

= 4 



' ^ 2 2 



350 keV 
15.1 keV 
320 keV 
50-1000 keV 
263 keV 
if keV 
13.4 keV 
2.6 MeV 
560 keV 
6.3 keV 



NR,hc 
NR,hc 
NR,hc 
NR,FF 

R 
R,IM 
NR,hc 
NR,FF 
NR,hc 
NR,hc 











Unfortunately, the simple theory says that Be (of which the beam pipe is made) 
will not bind to monopoles {S = 3/2, 7 < 0), but the strength of the magnetic 
field in the vicinity of a monopole will mostly likely disrupt the nucleus, ensuring 
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binding. Al and Pb° are okay. 

Estimates of binding energies are in the keV range and up. Simple estimations 
show that an energy of an eV is sufficient to bind a monopole to a nucleus with a f 
year lifetime, and then the monopolc-atom complex will remain permanently bound 
to the material lattice. (Evidently, it would not do to melt the material down.) 

A characteristic signal would be produced in a superconducting loop contained 
within a superconducting can by a magnetic monopole of strength g pulled through 
it. Note that if the shield were not present, the supercurrent induced in the loop of 
inductance L and radius R would be given by 

m-^-^(i^^^^=], (46) 

where z{t) is the vertical position of the monopole relative to the position of the 
center of the loop. The theory including the shield correction can be verified with 
a pseudopole. 

Background effects are enormous. All nonmagnetic but conducting samples 
possess: 



Permanent magnetic dipole moments (in the absence of boundaries) : 

m = _ ^ ^ (47) 

^' LC [i?2 + z(t)2]3/2 y'^') 

Induced magnetization: (a is the radius of the superconducting cylinder) 

V I f , , T , , dB, , f f z' a\ , ^ 

m--,-Jidr)rMr)^i^')j,H (48) 

where v is the velocity with which the sample is pulled through the detector, 
a is the conductivity of the sample, and 



OO 



Ki[y)-h[v) 



hiya/R) 



(49a) 



i?3 

a/R oo. (49b) 



2 (i?2 + z'^)3/2 



Calibration and real data are shown in the Figures. The pseudopole data (Fig. ^) 
clearly shows that we could detect a Dirac pole. As one sees from Fig. ^ real 
samples have "large" dipole signals; what we are looking for is an asymptotic "step" 
indicating the presence of a magnetic charge. Steps seen are typically much smaller 
than that expected of a magnetic pole of Dirac strength. The histogram of steps 
is shown in Fig. ^j. For n = 1 the 90% confidence upper limit is 4.2 signal events 

"22% of naturally occurring Pb is ^gJPb, whichibfis spin 1/2 and 7 = 0.582, which is sufficient for 
binding in the lowest angular momentum state. t3 The other three stable isotopes have spin 0. 



Milton, Kalbfleisch, and Gamberg 



c 
o 

Q. 

m 

0) 



> 

0) 

tn 

o 

Q. 

0) 
DC 



0.8 - 



0.4 



A A 



0.0 |-p,j5a&a&A^' 



-30 



-20 



H h- 

-10 
Position 



10 

5 

-5 
-10 



10 

(cm) 



la 

Monopole 

■ Calib. Data 

Dipole 

* , , Dipole Data 



20 



30 



■ 2.3 poles 

- Step fit 

5.5 mV 


lb 

t 


---..;:■■"■/ 

— 1 1 ' 1 





20 



40 

Position 



60 
(cm) 



80 



100 



Fig. 3. "Pseudopole" curves, a) Comparison of theoretical monopole response to an experimental 
calibration and of a simple point dipole of one sample with that calculated from the theoretical 
response curve, b) The observed "step" for a pseudopole current, corresponding to 2.3 minimum 
Dirac poles, embedded in an Al sample. 



Table 2. Acceptances, upper cross section limits, and lower mass limits, as determined in this 
work (at 90% CL). 

Magnetic Charge \n\ = 1 |n| = 2 |n| = 3 \n\ = 6 

Sample Al Al Be Be 

An/47r acceptance 0.12 0.12 0.95 0.95 

Mass Acceptance 0.22 0.060 0.0018 0.11 

Number of Poles < 4.2 < 2.4 < 2.4 < 2.4 

Cross section limit 0.84 pb 1.9 pb 8.4 pb 0.17 pb 

Monopole Mass Limit > 263 GeV > 282 GeV > 284 GeV > 413 GeV 



for 8 events observed when 10 were expected. These 8 samples were remeasured 
and all fell within ±1.47 mV of n = (more than 1.28cr from n = ±1). For n = 2 
the 90% confidence upper limit is 2.4 signal events for zero events observed and 
zero expected. Then, by putting in angular and mass acceptances we can get cross 
section limits as shown in Table ^ These numbers reflect a new analysis, and so 
differ somewhat from our published results. Ell To obtain the mass limits, we use the 
model cross sections referred to above. 
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regions of the plot in this vertically expanded view. 
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